Cosmic recall and the scattering picture of Loop Quantum Cosmology 
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The global dynamics of a homogeneous universe in Loop Quantum Cosmology is viewed as a 
scattering process of its geometrodynamical equivalent. This picture is applied to build a flexible 
(easy to generalize) and not restricted just to exactly solvable models method of verifying the 
preservation of the semiclassicality through the bounce. The devised method is next applied to 
two simple examples: (i) the isotropic Friedman Robertson Walker universe, and {ii) the isotropic 
sector of the Bianchi I model. For both of them we show, that the dispersions in the logarithm of 
the volume \n{v) and scalar field momentum \n{p^) in the distant future and past are related via 
strong triangle inequalities. This implies in particular a strict preservation of the semiclassicality 
(in considered degrees of freedom) in both the cases (i) and (ii). Derived inequalities are general: 
valid for all the physical states within the considered models. 

PACS numbers: 04.60.Pp,04.60.Kz,98.80.Qc 



I. INTRODUCTION 

Loop Quantum Gravity [H, [3| and its symmetry re- 
duced analog, known as Loop Quantum Cosmology [1, Q 
have experienced over recent years a dynamical progress. 
In particular, an application of the latter to the stud- 
ies of a simplest (isotropic) models of an early Universe 
have shown that the quantum nature of the geometry 
qualitatively modifies the global picture of its evolution. 
Namely, the big bang singularity is dynamically resolved 
as it is replaced by a so called big bounce @ which con- 
nects the current (expanding) Universe with a contract- 
ing one preceding it. The results initially obtained nu- 
merically for the flat isotropic model with massless scalar 
field_^,[71] were shown to be general features of that model 
and next extended to more general matter fields 
■|12| and topologies as well as to less symmet- 

ric systems: anisotropic ^15i - i21i | and some classes of the 
inhomogeneous ones Also, although the theory orig- 
inates from the canonical framework, a connection with 
the spin foam 23J models was made through the studies 
of the path integral in LQC [l^l as well as the analysis 
of the cosmological models within the spin foam mod- 
els themselves . Another avenue of extensions is the 
perturbation theory around the cosmological solutions 
|26| . The elements of its mathematical structure of LQC, 
which was initially formulated in 27[, were investigated 
in detail fH^Hl • In addition to the studies performed on 




the genuine quantum level, there is a fast growing num- 
ber of works einploying a classical effective formulation of 
the dynamics [32l, l33| , which often provides qualitatively 
new predictions [34,~36i] . The methods of LQC were also 
applied, with various levels of rigor, outside of the cosmo- 
logical setting, in particular in description of black hole 
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solutions [371, |38| and spherically symmetric spacetimes 
[39I lioj . There exist also studies of different prescriptions 
within the polymeric quantization ^4T1 - |43| as well as of 
the connection between LQC and the noncommutative 
geometry [iS]. The effects predicted by LQC were also 
applied for the regularization of the cosmological models 
not originating from the polymer quantization ^Bj. 

The prediction of the bounce and the existence of the 
branch of the universe evolution preceding it have raised 
an interesting question: provided that the expanding 
post-bounce branch is semiclassical, what can we deduce 
about the pre-bounce one? Does it have to be necessarily 
semiclassical or can it be completely dispersed and not 
possible to describe by any classical metric? The prelim- 
inary studies performed in context of the simplification 
of LQC and presented in [i^ seemed to favor the latter 
possibility. However more detailed analysis performed 
in the framework of the so called solvable LQC [1] have 
shown, that for the states satisfying quite mild semiclas- 
sicality assumptions for one of the branches the possible 
growth of the dispersion through the bounce is severely 
limited In consequence, within considered class of the 
states the semiclassicality at the distant future implied 
the semiclassicality in the distant past [47| . 

The latter results, although firm, strongly rely on the 
analytic solvability of the studied model, thus are difficult 
to extend to the original formulation of it [7] as well as 
to more general settings, for which one could not a priori 
exclude the loss of semiclassicality through the bounce 
p8| . Finding a definitive answer to the question posed 
above requires construction of the method, which is suf- 
ficiently flexible to be adaptable to the situations, where 
the analytical studies fail. We introduce such method 
in this article, next applying it to two simple examples 
of flat models with massless scalar field as the sole mat- 
ter content: [i) an isotropic Friedman-Robertson- Walker 
universe, and {ii) an isotropic sector of a Bianchi I model 
quantized as specified in [20| • In both cases our technique 
allows to derive certain (strong) triangle inequalities in- 
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volving the dispersions of (the logarithms of) the total 
volume and the scalar field momentum. The inequalities 
are general, valid for all the physical states admitted by 
the model. In the case (z) they imply an exact preserva- 
tion of the semiclassicality, whereas in (m) an analogous 
result is only partial, as it does not involve all the degrees 
of freedom describing the anisotropic system. 

The treatment we introduce here is based on the ob- 
servation, that the structure of the evolution operator 
in LQC implies that for each physical state (in cer- 
tain sense) there exists the geometrodynamical (Wheeler- 
DeWitt) one, which is its large scale limit in either strict 
0j [HI or approximate [lol . [l3| sense. Furthermore in a 
large class of the models the structure of the physical 
Hilbert space of each geometrodynamical theory describ- 
ing the limit admits a decomposition onto equivalents of 
Klein- Gordon plane waves either incoming from or out- 
going to an infinite volume. Since the WDW limits of 
LQC states are formed out of "standing waves" coupling 
both the incoming and outgoing ones, one can perceive 
the LQC evolution as a scattering process, which trans- 
forms the WDW incoming state (contracting universe) 
onto the outgoing one (expanding). Known properties 
of the limit allow to explicitly determine the scattering 
matrix and thus to relate the properties of the incoming 
and outgoing states. The procedure is explained in detail 
in Section HVl through an application to the cases (z) and 
{ii) listed above. 

The paper is organized as follows: First, in Sec. HI] 
we briefly introduce the main aspects of the LQC frame- 
work used to characterize considered models, as well as its 
WDW analog. Next, in Sec. IIIII we analyze in detail the 
exact WDW limit of LQC states, which is next used in 
Sec.|lV]to construct the scattering picture mentioned in 
the previous paragraph. That picture is then employed 
in Sec. |V] to the analysis of the dispersions of the in- 
coming/outgoing WDW states, which allows in turn to 
derive the triangle inequalities relating them. We also 
show there fSec. lVBj) . that the dispersions of the outgo- 
ing/incoming components of the WDW limit equal the 
dispersions of a genuine LQC state in the asymptotic fu- 
ture and past, thus extending the inequalities to these 
quantities. We conclude with Section IVII with the dis- 
cussion of the main results as well as the possibilities 
of their extensions to more general settings. In order to 
make the presentation of the studies clearer to the reader, 
some details of the mathematical studies as well as the 
details of the numerical analysis used in the article are 
moved to Appendix [X] and |B] respectively. 



II. THE FRAMEWORK OF LQC AND ITS 
WDW ANALOG 

In this section we briefly sketch the specification, quan- 
tization program and relevant properties of the models 
we are going to study. We discuss only those elements of 
the theory, which are relevant for our analysis. For the 



detailed description of the quantization and the proper- 
ties of the systems the reader is referred to 0, [131 and 

Since the considered systems are constrained ones, 
they are quantized via Dirac program, which consists 
of the following steps: kinematical quantization ignoring 
the constraints, promoting the constraints to quantum 
operators and finding the physical Hilbert space built out 
of the states annihilated by the quantum constraint. In 
accordance to this procedure this section has the follow- 
ing structure: First we introduce the classical models, 
LQC kinematics and quantization of the constraint for 
the FRW and Bianchi I models separately in Sec. Ill Al 
and |TTB] respectively. Next we describe the structure of 
the physical Hilbert space and relevant observables for 
LQC theory (in Sec. IIIC|) . as well as its WDW analog 
(SeclHDl). 

A. Isotropic flat FRW universe 

1. Classical theory 

On the classical level spacetimes of this class admit 
a (parametrized by a time t) foliation by homogeneous 
surfaces A4 — x M., where E is topologically R"^. Their 
metric tensor is 

g = -N^dt'' + a^{t)''q, (2.1) 

where is a lapse function, a is a scale factor (or equiv- 
alently a size of certain selected region V, see the discus- 
sion below (|2.3|) ) and °q is a fiducial Cartesian metric. 
To describe the system further we apply the canonical 
formalism, expressing the geometry in terms of Ashtekar 
variables: connections and triads, selecting the gauge fix- 
ing in which they can be expressed in terms of the real 
connection and triad coefficients c,p 

where ce^ / °uj\ is an orthonormal triad/cotriad corre- 
sponding to the fiducial metric °q and Vo is the fidu- 
cial volume of V. The variables c,p are canonically con- 
jugated with {c,p} — 87rG7/3 (where 7 is a Barbero- 
Immirzi parameter, which value has been set following 
[49} ) and are global degrees of freedom of the geometry. 
In particular p = a? . 

Within selected gauge all the constraints except the 
Hamiltonian one are automatically satisfied. The remain- 
ing constraint takes the form C — N{C^^ + C^) where 

C^r^-\l d^xe^.kC-'E^^E'^F^,, (2.3) 
7 Jv 

with e := det(£')| and F being a curvature of a con- 
nection A. To deal with the noncompactness of E the 
integration of a Hamiltonian density was performed only 
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over a chosen cubic cell V constant in comoving coor- 
dinates, which is an equivalent to the infrared cut-off. 
Despite this, the physical predictions are invariant with 
respect to the choice of the cell fl). 

The remaining matter part of the constraint equals 

C4. = 8nGpl/p'/^ (2.4) 

where (j) and are, respectively, the value of the scalar 
field and its conjugate momentum with {(j),p^} = 1. 



2. Loop quantization 

The classical system specified in Sec. Ill A II is next 
quantized via methods of Loop Quantum Gravity. In 
particular, the Dirac program is employed to construct 
the physical Hilbert space. It consists of the following 
steps: 

• Quantization on the kinematical level (ignoring the 
constraints). Here, the as the basic objects instead of 
A^, we select the holonomies of along the straight 
lines and fluxes of Ef along the unit square 2-surfaces, 
which form a closed algebra. The direct implementation 
of the procedure used in LQG [2i] leads to the gravita- 
tional Hilbert space = L^(]R, d/iBohr), where R is 
the Bohr compactification of the real line. The basic 
operators are respectively the holonomies /i^"^^ along the 
edge of the length A and triad p (corresponding to the 
flux across the unit square). The basis of "Hgr is built 
of the eigenstates of p and parametrized by v such that 
p\v) = (27r7€p[-\/A)^^'^ sgn(?;)|wp/'^|u), where the parame- 
ter A is the so called area gap specified in the next point. 
The inner product on Hgr is given by 

The matter degrees of freedom are quantized via stan- 
dard methods, thus attaining the Schroedinger-like rep- 
resentation. In consequence the full kinematical Hilbert 
space takes the form 

Hkin = ® n^, i'(M, d0). (2.6) 

The basic operators on are the field value (f> and its 
momentum p^ = —ihd^. 

• Promoting the constraint to quantum operator. For 
that all the geometric components in (|2.3|) and ()2.4p have 
to be expressed first in terms of the holonomies and 
fluxes, which is essentially done via methods specified 
in [5^. The field strength is in particular repre- 
sented via holonomies along the closed square loop. The 
requirement for our theory to mimic the properties of 
LQG and the discreteness of the area operator Ar there 
forces us to fix the physical area A of this loop as the 1st 
nonzero eigenvalue of Ar, which is the unique physically 
consistent choice for that technique |35l] . 



Presently in the literature there exists several pre- 
scriptions of constructing the quantum Hamiltonian con- 
straint differing by fine details, like the lapse, the factor 
ordering and the symmetrization. Here we focus on three 
of them, defined in Q, [1| and [H, HII and denoted, re- 
spectively, as the APS, sLQC and MMO prescription. In 
all these cases the resulting operator can be brought to 
the form 

t®dl + Q®t, (2.7) 

where an action of the operator Q equals 

-[Q^]{v) = f+(v)i,(v - 4) - u^mv) 
+ fAv)^{v + ^), 

with the form of fo^± depending on the particular pre- 
scription used and given respectively by 

• APS: 

Uiv) = [B{v ± 4)]-i/(t; ± 2)[Biv)]-K (2.9a) 
Uv) ^ [Biv)]-'[U{v) + Uv)], (2.9b) 

where [52] 

f{v) ^ {3nG/8)\v\ I |w + 1| - |w - 1| I (2.10a) 
Biv) = (27/8)|w|||w + l|i/3_ |„_i|i/3|3 (2.10b) 

• sLQC: 

f±iv) = i37rG/4)^v{v±i){v ± 2), (2.11a) 
fo{v) = (3^G/2)z;2, (2.11b) 

• MMO: 

f±{v) - Cg{v ± 4)s±{v ± 2)g^{v ± 2)s±{v)g{v), 
fo{v) = Cg\v)[g^{v - 2)sl{v) + g^v + 2)4(«)], 

where 

g{v) = ||1 + l/?;|i/3 _ |i _ l/«|i/3|-i/2^ (2.13a) 

s±(u) = sgn(t;±2) +sgn(w), (2.13b) 

C^nG/U. (2.13c) 

The operator 6 is denoted, respectively, by 6aps, ®sLQC 
and 0MMO a-iid is well defined in all the listed prescrip- 
tions in particular for e = (see the detailed discussion 
in 131] for APS and [51] for MMO). 

• Building a physical space out of states annihilated 
by the constraint. Since the operator (|2.7p is essentially 
self-adjoint [11], this step can be realized via system- 
atic procedure of the group averaging [5^ [13] • On the 
other hand its form selects another natural way of find- 
ing the solutions 0, 113] which in this case is equivalent 
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to it, namely the reinterpretation of the constraint as the 
Klein-Gordon-like equation 



(2.14) 



defining the evolution of a free system along 0, that is 
the mapping between the the spaces of the "initial data" 
- restrictions of ^ to the surfaces of constant (h 



3 (/.K^ *(•,</>) e -Hg 



(2.15) 



The quite simple form of the operator Q allows to eas- 
ily define the physical Hilbert space "Hphy via its spectral 
decomposition. This step, as well as the notion of evolu- 
tion will be described in more detail in Sec. Ill CI 

Before going to it let us note, that the structure of Q 
and (|2.7p provides the natural division of the domain of 
V onto the subsets (the lattices) 



^{e + An; n£ Z}, e e [0,4[ 



(2.16) 



preserved by the action of Q. This division is naturally 
transferred to the splitting of "Hphy onto the superselec- 
tion sectors. In consequence it is enough to fix particular 
value of e and work just with the restriction of the do- 
main of 6 to functions supported on only. For the 
clarity we will consider just the sector e = 0, however 
the presented treatment and its results generalize easily 
(at the qualitative level) to all the sectors. 

Further simplification comes from the fact, that the 
considered system does not admit parity violating in- 
teractions. In consequence the triad orientation reflec- 
tion V t-^ —V being the large symmetry provides another 
natural division onto superselection sectors, namely the 
spaces of symmetric and antisymmetric states. For the 
selected sector e = this particular choice allows to 
further restrict the support of the functions to £q 



B. Flat Bianchi I universe 

The first step in the generalization of the model pre- 
sented in previous Section is an extension to the flat 
Bianchi I model, describing the universe with the same 
matter content and topology, which however while being 
still homogeneous is not necessarily isotropic. Its (pre- 
liminary) analysis within LQC has been initiated in |15| . 
Later more detailed analysis of its kinematics and dy- 
namics was performed in 16 1 and [13] (see also [isl. Il9j 
for a vacuum case), although the quantization prescrip- 
tion used there is not applicable to the noncompact cases 
[361] . The first description valid also in noncompact sit- 
uation was constructed in [20], which we will follow in 
this article. In this section we briefiy introduce those 
elements of the framework, which are needed as a basis 
for our analysis. The treatment is in fact an extension 
of the one applied to the FRW universe. Therefore, for 
shortness, here we will focus just on these aspects of it, 
which differ from the description presented in Sec. Ill Al 
For the detailed description the reader is referred to j2^] . 



1. The classical model 

Classically the Bianchi I fiat spacetime admits the 
same foliation by homogeneous surfaces as the FRW one. 
The general form of the metric is 

g = -N^dt^ +al{t)dxi^ + al{t)dx2^ + al{t)dx3^ (2.17) 

where Xi are chosen (comoving) coordinates defining the 
fiducial Cartesian metric °q := da;i^-|-da;2^-|-da;3^ withor- 
thonormal (co)triad {°ujD oef . To construct the Hamil- 
tonian formalism one has to introduce again a fiducial 
cell, which here is described by three fiducial (i.e. with 
respect to the metric °q) lengths Li and fiducial volume 
Vo '■— L1L2L3. Here it is again possible to fix a gauge 
in which the Ashtekar connections and triads are rep- 
resented by three pairs of canonically conjugated coeffi- 
cients c\pi 
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and all the constraints except the Hamiltonian one are 
automatically solved. The Poisson brackets between co- 
efficients equal {d,pj} — SnGjS'j. 

Following [l^l we choose the lapse N = y/\piP2P3\- 
The Hamiltonian constraint has the same form as in the 
FRW case, and its components are given by (|2.3p and 
(with the basic variables given by (|2.18p instead of 
])), where in p := \piP2P3\^^^- 



2. Loop quantization 

To quantize the system we follow the program spec- 
ified in Sect. Ill A 21 in particular choosing the polymer 
representation for the geometry degrees of freedom, while 
keeping the Schroedinger one for the scalar field. 

In the geometry part the basic operators are 
holonomies along straight edges generated by oe^ and 
the fluxes across 2-dimensional rectangles spanned by 
oef . The gravitational kinematical Hilbert space consists 
of the product of three copies of Hgr of the FRW sys- 
tem, each corresponding to one direction of oe° Hgr — 
0^^;^ L^(K, d/ZBohr)- The basis of this space can be built 
out of eigenstates of the triad (or unit flux) operators 
Pi and parametrized by three real variables \i such that 
ft]Ai,A2,A3) = sgn(A,)(47r7v^4i)lA|]Ai,A2,A3). Al- 
ternatively, one of A^ can be replaced with the parameter 
V := 2A1A2A3. Here for that purpose we select A3, finally 
labeling the basis elements as ]Ai, \2tV). 

The space and the set of basic operators corre- 
sponding to the matter are the same ones as in the FRW 
case. The full kinematical Hilbert space is also a product 

■Hkin = Hgr ® H^. 

The quantum Hamiltonian constraint is constructed 
out of the classical one by, first reexpressing it in terms 
of the holonomies and fluxes, and next promoting these 
components to operators. In particular the field strength 
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^ab again represented via holonomies along closed rect- 
angular loops of the physical area equal to A. Unlike in 
isotropic case however fixing the loop area does not al- 
low to uniquely fix the fiducial lengths of its edges. This 
apparent "ambiguity" gave rise to several distinct pre- 
scriptions present currently in the literature (including 
the one of On the other hand the relation of the 

LQC degrees of freedom with the full LQG ones con- 
structed in [l^l allowed to fix the relation uniquely. The 
constraint resulting from this operation (defined on the 
dense domain in Hkin 'SS' '^(j?) is of the Klein-Gordon form 

(7=1(^9^ + 1(^631. (2.19) 

Similarly to FRW model we can restrict our interest 
to just the symmetric sector, that is those states 4', 
which satisfy *(Ai, A2, t", </)) = ^'(|Ai|, IA2I, 0). This 
allows to restrict the studies just to the positive octant 
Ai, A2, u > 0, on which an action o f 8b 1 is given by (quite 
complicated) Eqs. (3.35)-(3.37) of [20]. Its important fea- 
ture is, that, analogously to the isotropic one, it divides 
Hphy onto the superselection sectors built of the states 
supported just on the sets {(Ai, A2, f ); Ai, A2 G M, u e 
£+} with £+ := {v = e + 4n; n € N}, preserved by an 
action of Qbi- Therefore to extract the physics one can 
consider just one of those sectors. Here, for simplicity we 
choose e = 0. 

As we show in Appendix I A II the operator 6bi admits 
self-adjoint extensions. Knowing its action one can in 
principle find the physical Hilbert space(s) corresponding 
to the model by analyzing the spectral properties of the 
extensions. On the other hand there exists a well defined 
procedure of the averaging over anisotropics (defined in 
[20"] to build an embedding of the isotropic model in the 
homogeneous anisotropic one). In this article we will 
focus just on the space of averaged states Uphy and their 
physical properties. 

3. The isotropic sector 

Following [l^] and the ideas of [11] we consider a pro- 
jection P mapping from the dense domain in "Hgr of the 
Bianchi I model to the one the isotropic model as follows 

V'(Ai,A2,«) ^ ^(w) = [Pij]{v) ^(Ai,A2,«). 

Ai ,A2 

(2.20) 

Through the direct inspection one can check that there 
exists an operator Obi such that 

[P0Bi^](«) ^ObiIPV-K^). (2.21) 

An action of that operator equals exactly the one of 
©sLQC defined via ([TgllCTl . 

The consequence of the above observation is, that at 
least to some extent those of the aspects of the Bianchi 
I model, which are related exclusively to the behavior 
of the isotropic degrees of freedom, can be investigated 



via an isotropic model constructed via averaging over 
anisotropics and equivalent to the sLQC one described 
in Sect. Ill Al However, one has to be aware, that some 
of the physical states might in principle be in the kernel 
of the projection operator P. Thus, certain care needs 
to be taken, when relating the properties of the isotropic 
sector defined above with the full Bianchi I model. In 
particular it is not confirmed, whether the expectation 
values and dispersions of the total volume Bianchi I op- 
erators agree with the analogous quantities of the volume 
operators acting on the averaged states. This issue will 
require further studies. 

C. Physical Hilbert space, observables 

Known form of the Hamiltonian constraint (j2.7l) and in 
particular the evolution operator O (|2.8p allows to easily 
extract the Hilbert structure of the space of states an- 
nihilated by the constraint. The exact construction of 
Hphy is done via group averaging (see for the details) . 

To start with, we note, that the spectrum of is for 
all the cases considered here absolutely continuous [56| . 
nondegenerate and equals Sp(0) = R"*" U |0| (29j . In 
consequence one can build a base of Hgr [571 ] out of the 
eigenfunctions ek corresponding to nonnegative eigenval- 
ues 

[eek]{v)^uj^ik)ek{v), (2.22) 

(where uj{k) = y/T2nGk, fc > 0) and normalized such 
that (efc]efc') = 6{k' — k). In the superselection sectors e = 
the remaining freedom of global rotation is furthermore 
fixed by the requirement, that 

efc(u = 4)eR+. (2.23) 

Applying the simplest form of the group averaging pre- 
sented in [H and above spectral decomposition we arrive 
to the following representation of the elements of Hphy 

^-(^,0)=/" dk^{k)ek{v)e"^^''^'f', (2.24) 
Jr+ 

where e L'^(R+,dk) is the spectral profile of The 
physical inner product is given by 

(^f]$) = / dfc*(A:)$(/c). (2.25) 

Jr+ 

Since here we are dealing with the constrained system, 
there is no natural notion of time and evolution. It can 
be provided via the unitary mapping (j2.15p . An alterna- 
tive way to define an evolution is the construction of the 
family of partial observables [11], parametrized by one of 
the dynamical variables and of the elements related via 
unitary transformation. Here we construct family In]?)]^ 
16;] interpreted as In ]u] at given "time" 0. The systematic 
way of constructing such observables is presented in [3(| . 
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For the models considered in the article the expectation 
values and dispersions of In \v\^ for the physical state ^ 
equal the analogous quantities of the kinematical observ- 
able In j-OI acting on the initial data "0^ :— ^(-,0) S Hkin 

[ln(i;)^„*](w) = e^^f-^--^") \n{v)^{v, (j^o). (2.26) 

For completeness we introduce one more observable, cor- 
responding to the constant of motion ln(a;) - an operator 
In(ii) acting as follows 

[ln(w/%/G)VE'](A:) = ln(a;(fc)/v^)*(fc). (2.27) 

This operator will be useful later in the paper as (it will 
be shown that) its dispersion bounds the growth of the 
spread in In jSj^. 



D. Wheeler-DeWitt analog 

The systems studied in this article can be also quan- 
tized via methods of the geometrodynamics. Indeed, the 
geometric component (j2.3l) of the Hamiltonian constraint 
can be expressed entirely in terms of the coefficients (c, p) 
defined in 

C,,=-^^'Vp, (2-28) 

and the entire system can be treated just as an abstract 
one of the phase space coordinatized by {c,p,(p,p^) and 
quantized via standard methods of quantum mechanics. 
As the result the kinematical Hilbert space takes the form 
2ikin = Kgr®-^^(I^. df/*), where Hg, L'^{R, dv) with the 
inner product between V', X € Hgr 

{±\X) = I dv±{v)x{v). (2.29) 

The quantum Hamiltonian constraint can be expressed 
as a differential analog of (12. 7|) 

To arrive to above equation we selected the factor order- 
ing in (j2.28p consistent with the one of (j2.3p (see j^l for 
the detailed explanation). 

The physical Hilbert space can be again constructed 
via group averaging and it is an almost complete analog 
of the one of LQC models, with just slight differences 
being a consequences of the two-fold degeneracy of the 
eigenspaces of the operator 0. Here the orthonormal 
basis of (the symmetric states on) Tigr consists of the 
functions 

e^{v) ^ ^=e^^'^^\^\ fceM, (2.31) 



and the physical states (positive frequency solutions to 
([230]) ) have the form 

3r(„,0)= j dfct(fc)efe(w)e^"'('=)'^, (2.32) 

where * e L'^iR.dk) and uj{k) = Vl27rG|fc|. The inner 
product has the same form as ()2.25p but now k runs the 
entire R. 

To characterize the states and define a physical evo- 
lution we use the observables ln(w) and In |t}|0 given, re- 
spectively, by full analogs of ^[TTfl and (P^ . The latter 
ones can be expressed as quite simple differential opera- 
tors acting directly on Kphy 

In Iwl^i = -ie*"('=)'^9fce-*"('=)*£ 

= [-idk - {dMk))4>i]±. ^^'^^^ 

This fact will be very useful in the following sections, 
where we will use it to derive the relation between the 
dispersion of the components of the WDW limits of the 
LQC states. 



III. WDW LIMIT OF AN LQC STATE 

The comparison of the forms of the operators Q (|2.8p 
and 6 specified via (|2.30p shows that under certain con- 
ditions (slowly changing functions) one of the operators 
can be approximated by the other. Therefore one may 
expect, that the solutions to (12.71) converge in certain re- 
gions to some solutions of (I2.30p . Indeed, it was shown 
via numerical methods in [7] that the eigenfunctions 
converge to certain combinations of e^, and e_f^. Further- 
more the analytic properties (reality) of imply that this 
limit has the form of a "standing wave" that is it is com- 
posed equally of incoming {k > 0) and outgoing (fc < 0) 
plane waves (|2.31|) . More precisely [s^ 

e,{v) = r(A:)(e-('=)e,(^>) + e-»('=)e_,(«)) 

+ Oi\e,iv)\ik/vr) (3.1) 
=:l^{v) + 0{\e^{v)\{k/vf), 

where r{k) G can be determined analytically via the 
relations between norms of the LQC and WDW states 
(see Appendix I A 2|) and a{k) G 5^ is a phase shift. 

In this section we analyze the WDW limits V^^^ of the 
LQC eigenfunctions ek in detail. First, in Sec. IIII Al we 
provide an analytic proof of the convergence for all the 
forms of Q considered in this paper, as well as recall 
the arguments allowing to determine the structure of the 
limit. Second, in Sec. IIIIBI we perform a detailed an- 
alytical and numerical analysis of the phase shifts a{k) 
defined in (13.11) . The properties of these shifts are the 
critical components allowing to arrive to the triangle in- 
equalities relating the dispersions of the physical state 
and being the main result of this paper. 
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A. The convergence of the bases 

In order to explicitly show the convergence p.l[) we 
compare the eigenfunctions - solutions to (j2.22[) . with 
the solutions (IOT|) to the WDW analog of (g^H)- To 
start with, we note, that the Eq. (|2.22p is a 2nd order dif- 
ference equation, however, due to the specific properties 
of 6 characteristic for each of the prescriptions consid- 
ered here, the whole solution is determined just by the 
single value {v — 4) (see 0, H, [Ml for the details on 
respective prescriptions). 

To analyze the solution it is more convenient to rewrite 
that equation in the 1st order form ^oj. For that we 
introduce the vector notation, defining 



efc(w-4) 



(3.2) 



Using it we can rewrite (|2.22p in the following form 

+ A{v)ekiv), (3.3) 

where the matrix A can be expressed (with use of the 
notation introduced in Eq. (|2.8p ) as 



Aiv) 



/+(") f+{v) 

1 



(3.4) 



To relate Cfc with e^^ we note that the value of at 
each pair of consecutive points w, w -f 4 can be encoded 
as values of the (specific for the chosen pair of points) 
combination of e^f. , that is 



ek{v) = B{v)xk{v), 



(3.5) 



where the transformation matrix B is defined as follows 

(3.6) 



f efe(« + 4) efe(« + 4) 



Having at hand the objects defined above we can rewrite 
the equation (13. 3|) as the iterative equation for the vectors 
of coefficients Xk 



Xk{v + 4) = B-\v)A{v)B{v - A)xk{v) 
=: Miv)xk{v). 



(3.7) 



The exact coefficients of the matrix M{v) can be calcu- 
lated explicitly for each of the prescriptions specified in 
Sect. Ill A 21 An important feature (found by direct in- 
spection) of it is, that in all three cases they have the 
following asymptotic behavior 



M{v) = l + Oiv-^), 



(3.8) 



where 0{v~") denotes the matrix, all the coefficients of 
which behave like 0{v~"). That convergence implies (via 
application of the methods presented in ^ Sec. 4) the 
existence of the limit 



lim Xfc(4") = Xk, 



(3.9) 



as well as it confirms the rate of convergence specified 
in p.ip . This limit can be expressed in terms of the 
coefficients introduced in p.ip in the following form 



Xk = r{k) ( 



(3.10) 



which is a consequence of the observation, that all the 
coefficients fo,±{v) in (I2.22p are real, so (by (j2.23p ) is e^. 

The scaling factor r(fc) can be easily determined from 
the relation between norms in LQC and WDW theory 
discussed in Appendix lA 21 and equals 



(3.11) 



r(fc) = 2. 



The behavior of the phase shift function a{k) is however 
much less trivial and requires detailed studies. 



B. The phase shifts 

To extract the properties of the phase shifts a{k) de- 
fined in (I3.10p we combine the analytical and numerical 
methods. We focus on the behavior of the derivative 
a' := dka, as it is exactly the quantity which will be 
relevant in the further studies. First we derive analyti- 
cally the asymptotic behavior of a{k) for fc — > oo for the 
sLQC prescription (Sec. IIII B ip . The analytical results 
are then strengthened and generalized to other prescrip- 
tion by means of the numerical methods described in Ap- 
pendix |B] The results of that analysis are presented in 
Sec. HUB 21 



1. Asymptotics in sLQC 

Among the prescriptions considered here the sLQC one 
is somehow distinguished, as the bases Ck expressed as 
functions of an appropriately defined canonical momen- 
tum 6 of w have a simple analytical form [8] . This allows 
for a quite high level of control over their properties, a 
fact which we will exploit below. To start with, let us fix 
the definition of b, choosing it to equal 



boclpl" 



(3.12) 



where c, p are given by (|2.2p and the proportionality con- 
stant bo is fixed via the righthand side equality. This 
quantity can be next used as a configuration variable 
on the quantum level. A particularly convenient choice 
of the representation of the quantum states using this 
variable is provided by the following transformation op- 
erators chosen respectively for WDW and LQC {J-) 
framework 



dv \v\ 2 e" 'ip{v), 



(3.13a) 
(3.13b) 
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where the part of ip supported on < is determined 
by the symmetry requirement. The form of these trans- 
formations implies, that the domain of b is the entire R 
in the case of WDW and the circle of the radius 1 /2 in 
LQC. 

The evolution operators 8 and 6sLQC transformed via 
p.l3|) take the form 

= -12TTG[bdb]^ , (3.14a) 
©sLQC - -12^G[sin(6)4]2, (3.14b) 

and their (symmetric) eigenfunctions corresponding to 
the eigenvalue = 127rGfc^ (and in case of sLQC cor- 
responding to the sector e — 0) are combinations of the 
orthonormal basis elements 



=Fifeln|b/2| 



efe(6) =iV(fc)e-*'='"l''/2|, 

ek{b) = 7V(fc)cos(-fcln(tan|6/2|)), 



(3.15a) 
(3.15b) 



where N.{k),N{k) are the normalization factors deter- 
mined by the physical inner product [8] and (I3.15b[) is 
written in the chart b £ [— 7r/2, 7r/2]. The — sign in the 
exponents comes from the comparison of the spectrum 
and bases of the operator ^/Q in v and b representation 
(see for example |19j). 

To retrieve the large v behavior of the functions p.l5p 
one needs to perform the transform inverse to (j3.13l) . 
Since the large v correspond to high frequencies in 6, 
the particular form of the functions implies, that only 
the domain near 5 = will give the relevant contribution 
to the asymptotics in u oo. 

In order to be able to compare the functions e^, (b) and 
ek{b) one first needs to to deal with the fact, that the 
inverse transform of (j3.15bl) , involves an integration over 
the domain [— 7r/2,7r/2] whereas for p.l5al) one need to 
perform an integration over R. 

To do so, let us first consider on [~tt/2, tt/2] a function 
^(6)e;,(6), where ^(6) is some smooth function with sup- 
port in (— 7r/2, 7r/2) and equal to 1 in some neighborhood 
of 0. Regarding this function as defined on the entire R 
we note that the difference 



gk{b) -.^ e,{b) - mekib) 



(3.16) 



is a smooth function with appropriate behavior at infinity 
(6l| . Hence its Fourier transform is of the order 0{v^^) 
for any iV e N. 

On the other hand, for ^ (6)ej, [b] considered as a func- 
tion on a circle, the difference 



fk{b) 



-j:fcln(tan|f)/2|) 



me 



-ik\n \b/2\ 



(3.17) 



is of the class Cq , thus by Lebesgue-Riemann lemma its 
transform J-~^fk is of the order o(|ej,(u)|). In conse- 
quence the function 

e'kib) = m cos(fcln \b/2\) + (Mb) + f-k{b))/2, (3.18) 

supported on R has the same WDW limit as e^. Fur- 
thermore the components of this limit proportional to 



e^k correspond to the respective components e 



Bringing together this two observations we see, that 
to find the desired phase shifts one just needs to find 
the transform of the functions ej.(6). As they are the 
eigenfunctions of 8, they are proportional to e^iy) 



(3.19) 



and the factor of proportionality N{k) equals the trans- 
form of -v/27reT*''''"'''/^' at the points t; = ±1. Selecting 
for the component proportional to 6^,(6) the point v = —\ 
we get 



N{k) = V2^ / d6e-^('='"l^/2|-fc/2) 



(3.20) 



where to arrive to the latter equality we introduced the 
change of variables b = 2ky. The last integral can be 
computed in the regime /c — >■ (X) by a stationary phase 
method (see Appendix lA 31 for the proof of the applica- 
bility of the method). The result is 



N(k) « ^/8^fce^''='"('^)^/^2^ 



yo_^-ik[\n\y^\~yo) 

Vk 



Vo = l 



(3.21) 

Analogously, one can calculate N{k) for the component 
proportional to e_j.(6) by selecting the point v = 1. 
These two results allow us to extract the phase shift a(fc), 
which equals 



a{k) = -k(\n{k) - 1) - |7r + o(fc"). 



(3.22) 



Via the same method one can compute the derivative 
dkN =: N'. 



(3.23) 



N' = V2^ d6[iln|V2|]e-'('=i"l''/2|-''/2) 
Jr 

= y8^fce-*''-'"('=) J dy [-iln|y|]e-*'=^('"l''l-^) 

-\/8^ifcln(fc)e-*'= / dye-^'^C'^l^l-^) 
Jr 

w -i\n{k)N{k), 
where the last estimate follows from the fact, that de- 



composing N{k) =: A{k)e 



(where A{k) G R+) one 



can express its derivative as 

N\k) = ia'{k)A{k)e'-°'^^'^ + A{k)' e' 

« -iV8^\/^ln(fc)\/fce-'('='"('=)-'=) (3.24) 
+ 0{k-^'^\n{k)). 

In consequence the phase shift derivative equals 

a'(fc) = -ln(fc) +0(fc-Mn(fc)). (3.25) 
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2. Numerical generalization 

In the case of the remaining two prescriptions repeat- 
ing the analytical calculations preformed for the sLQC 
one is not possible, as the eigenfunctions of 8 do not have 
manageable analytic form in either of or 6 representa- 
tions. We note however, that between the prescriptions 
the operators & differ just by a compact operator. Thus, 
it is expected that the asymptotic behavior of both a{k) 
and a'{k) corresponding to them is again given by (I3.22p 
and (|3.25[) up to the rest terms decaying with k. We ver- 
ify this expectation for a'{k), using the numerical meth- 
ods which are described in detail in Appendix IbI Those 
methods allow to determine the values of a' in quite wide 
range of k as well as to verify its asymptotic behavior 
(within the limitations of applied numerics). The results 
for different prescriptions are presented in Fig. [TJ Al- 
though the exact form of a' depends on the prescription, 
especially for small fc, one can observe the following fea- 
tures common for all of the prescriptions considered in 
this article: 

(i) For large k the derivative a' converges to the limit 
specified in p.25p with the rate 

a'(fc) = -ln(fc) + 0(fc-2), (3.26) 

(ii) The scaled 2nd order derivative of a is bounded 

\[kdla]ik)\ < 1 (3.27) 

for every value of k. 

These properties will be crucial for building the relation 
between the dispersion growth through the bounce. 



IV. THE SCATTERING PICTURE 

It was shown in Sec. Illll that the basis functions span- 
ning the LQC physical Hilbert space admit certain WDW 
limits. Given that one can define a WDW limit of any 
physical state by replacing the basis functions in (j2.24|) 
with the limits ipj^ defined via ()3.1|) . This operation de- 
fines a relation between the LQC physical Hilbert space 
and the WDW one, which in terms of the spectral profiles 
can be written as follows 

^(|fc|) ^ i(fc) = 2e"sn(fe)a(|fc|) sgn{k)^{\k\), (4.1) 

where k spans the entire real line. 

That limit consists of two components, the incoming 
wave packet (corresponding to fc > 0) and the outgoing 
one {k < 0). On the physical level they represent the 
universe which is, respectively, contracting to big crunch 
singularity and expanding from the big bang one. The 
entire LQC dynamics can be thus seen as the specific kind 
of "transition" between the contracting WDW universe 



(represented by \^)- ) to the expanding one (denoted as 

^ l3^)out = Pli^)i„: (4.2) 

In consequence, looking at the evolution "from the infin- 
ity" (in the configuration space or in cosmic time) one 
can interpret the evolution as the process of scattering of 
the contracting geometrodynamical universe. The form 
of the limit (|4.ip immediately allows to write down the 
scattering matrix 

p{k,k') = (ejple,,) = e-^sn(fe'M|fe'|)^(fc^fc.)^ (43) 

which form encodes in particular the fact, that the con- 
tracting universe totally "reflects" into the expanding one 

±{k) ^ e2"Sn(fc')a(|fe'|)^(_^>)^ (^4 4^ 

This picture allows to address in a quite natural and intu- 
itive way the questions regarding the relation of the prop- 
erties of the pre bounce branch (universe in the asymp- 
totic past) and the post bounce one (asymptotic future). 
In particular, we will apply it in the next section to de- 
termine how much the bouncing universe can disperse in 
the distant future of the bounce in comparison to the 
initial spread in its distant past. 

When considering the above picture one has to remem- 
ber that, although the LQC basis functions converge to 
the combinations of the WDW ones, this is not neces- 
sarily the truth for the general physical states, as the 
convergence of the bases is not uniform with respect to 
k. Nonetheless, once the attention is restricted just to 
the states localized with respect to the observable k de- 
fined analogously to (j2.27p (that is of the finite dispersion 
in k) the uniformity is restored and the WDW limit is 
defined in the precise sense. This fact is used for example 
in Sec. IV Al where the expectation values and dispersions 
of the LQC states are related with the ones of its WDW 
limit. 

The scattering picture can obviously be constructed in 
the context of any LQC model in which the basis func- 
tions converge explicitly to the WDW ones, like for ex- 
ample the models with the positive cosmological constant 
|ll| or the Bianchi I ones [1^. The applicability is how- 
ever not restricted just to such systems. In particular the 
models featuring the classical recoUapse, like [lol. [l3|, in 
LQC admit a quasi-periodic evolution. For those models 
it is also possible to build a correspondence between the 
LQC and WDW states, since the basis elements of "Hphy 
converge to their WDW analogs also there. The new dif- 
ficulty in these cases is the fact that, as the spectra of 
the LQC evolution operators are discrete while the WDW 
ones are continuous, the direct analog of the transforma- 
tion (|4.ip leads to the WDW states of the zero norm. 
This problem can be circumvented by introducing an ap- 
propriate interpolation of the discrete spectral profile of 
an LQC state. The WDW state constructed this way 
represents a single epoch (between the bounces) of the 
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evolution of the LQC one. However, since such interpo- 
lations are not defined uniquely, there is no direct 1 — 1 
correspondence between the loop states and their "lim- 
its". One can however choose the interpolations which 
reproduce the expectation values and the dispersions of 
the relevant observables at least approximately. This way 
it is possible build the WDW states well mimicking one 
epoch of the loop universe evolution. It happens however 
at the cost of relaxing the relations between the physical 
parameters corresponding to them to approximate ones, 
without explicit convergence of their values. The reason 
for that is two-fold: 

(i) nontrivial corrections due to interpolation of the dis- 
crete spectral profile, and 

(ii) the fact that the basis elements of the LQC physi- 
cal Hilbert space converge to their (rescaled) WDW 
analogs only asymptotically, thus obviously beyond 
the classical recollapse point. 

Despite this, such relations can be still quite useful. In 
particular this method is well suited to address the ques- 
tion, how the parameters (for example dispersions) can 
change between the epochs. In particular it can be used 
to investigate the issue of the spontaneous coherence of 
the LQC state, that is to address the question whether, 
given an initial date describing the state which is not 
semiclassical, the state will admit in the future evolution 
the semiclassical epoch. 

V. THE DISPERSION ANALYSIS 

This section is dedicated to the main goal of this pa- 
per: finding the precise relation between the dispersions 
of the physical LQC state representing the universe in the 
distant future (post bounce) and past (pre bounce). The 
studies are divided onto two steps. First, in Sec. IV Al we 
apply the scattering picture to relate the dispersions of 
incoming and outgoing asymptotic WDW states. Found 
relation is next translated in Sec. IV Bl to obtain the re- 
lation between the dispersions of a genuine LQC state in 
the asymptotic past and future. 

A. Dispersions of the WDW limit components 

Given the WDW limit * (defined by dH])) of the LQC 
state described by the spectral profile let us define its 
decomposition onto the incoming and outgoing 
components such that the spectral profiles corresponding 
to them equal 

±^{k) ■.= 9{±k)±{k), (5.1) 

where 6* is a Heaviside step function. Denote the sub- 
spaces of 2iphy formed by these components as 2iphy re- 
spectively. On each of these components one can consider 



an action of the observables ln|t;|0 defined by (|2.33p . 
Their expectation values and dispersions equal respec- 
tively 

vi^) := {\n\v\^)± = ±a0||3^±|| + Hdk)±, (5.2a) 
a± := (A In |i;|^)± = {A{~zdk))±. (5.2b) 

where a :— vl27rG and for any observable O we define 
(6)±:=(*±|6|*±). 

The main question we would like to address here is 
whether there exists the relation between cr_ and a+ and 
what is its form. The answer to the former is certainly 
true as the transformation (|4.4I) unitarily maps 
in the following way 

±-{k) - [C/i+](fc) = e^2a(|fe|)|,+ (_^)^ 

+ _ (5.3) 

U ■ Kphy ^ 2iphy' 

thus the expectation values and dispersions in (|5.2p are 
related as follows 

{^tdk)- = {U-^Hdk]U) + , (5.4a) 

{A[-tdk])- = {AU-^[~idk]U) + , (5.4b) 

where 

U^^[-idk]U = -idk -2a'l. (5.5) 

Combining together (|5.2bp . (jOl) . ()5.4bp . ([53]) and ap- 
plying very general bound on the dispersion of the sum 
of operators (|A24[) we obtain the following inequality 

cr_ < (7+ + 2(Aa'l) + . (5.6) 

To write it down in the useful form we have to express the 
quantity 2(Aa'l)+ in terms of dispersions of observables 
commonly used to characterize the physical properties of 
the state. For that we exploit the properties of the func- 
tion a' found in Sect. IIIIBI Namely, by the definition of 
the dispersion we can bound the term under considera- 
tion via 

{Aa'l)l = {{a'' - {a') + f 1)+ < {{a'' - a'*f 1) + , 

(5.7) 

which is true for any value of a' . Here we choose it to 
be 

a'* = a'(cxp(A*)), A* := (ln(A:)) + . (5.8) 

Upon that choice, applying the property (|3.27|) of a we 
can bound the left-hand side of (|5.7p as follows 

{Aa't)\ < ((In(fc) - = (Aln(fc))^. (5.9) 

Finally, knowing the relation uj{k) we can express the 
right-hand side of (|5.9p via the dispersion of the WDW 
analog of the observable (|2.27l) . which corresponds just 
to a logarithmic scalar field momentum ln(p0/6), where 
b := hVC. The resuh is 

<a+ + 2a^, (5.10) 
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where 

a. = (Aln(p^)/6)+ = (*|Aln(p^)/6|v|/). (5.11) 

The righthand equaUty is here a direct consequence of 
the form of operator (|2.27p (multiphcation operator in 
k) and the transformation (14.11) . 

One can easily see, that the role of cr_ and (T+ can be 
exchanged. The only modification induced by this oper- 
ation will be the exchange U O U~^. In consequence, 
(|5.10p is supplemented by the inequality 

a+ < o-_ +2(7^, (5.12) 

thus both these inequalities (|5.10p and (|5.10l) can be un- 
derstood as the triangle inequalities. 

Note that to arrive to above inequalities we have not 
assumed any semiclassicality conditions in any epoch of 
the state evolution, neither we required the state to be 
peaked about any trajectory. The relations hold for every 
element of "Hphy 

It is also worth noting, that although (|5.10p . (|5.10p are 
formulated in terms of the dispersions of the logarithmic 
observables Injwl^ and ln{p^/b), for the states semiclas- 
sical (sharply peaked) in any epoch of the evolution (pre 
or post bounce) these quantities can be approximated via 
analogous "linear" ones: juj^, 

(In \v\^) « ln(|5U), (HP^/b)) « H{p^)/b), 

(Aln|^U)«i^, (ln(p,/fe)) « 

In consequence the inequalities (|5.10p . (|5.12p can be re- 
formulated in terms of them at least on the semi-heuristic 
level (or in precise sense under certain additional assump- 
tions imposed on the form of the state). This is not 
however the aim of the article, as we intended to find 
a relation which is maximally general while remaining 
precise. 

At this point one has to be aware of an important 
subtlety related to the description method applied hare. 
Namely, the considered observables are the geometro- 
dynamical observables acting on the asymptotic states 
(wave packets) , not the exact LQC observables acting on 
the LQC states. Therefore one may in principle worry, 
that due to some wild behavior of the LQC basis func- 
tions near the bounce point there might be some residual 
contributions to the results (expectation values, disper- 
sions) of the scattering picture essentially invalidating 
found inequalities, once applied to exact LQC observ- 
ables. 

On the other hand, the studies of Q performed for "lin- 
ear" observables \v\^ show, that at least for the sLQC pre- 
scription the LQC dispersions indeed approach the ones 
of WDW limits. This suggests that the problematic cor- 
rections mentioned above are not sufficient to distort the 
main results. Indeed, one can confirm this expectations 
in quite general setting using the relation of the norms of 
LQC state and its WDW limit derived in Appendix lA 21 



B. The relation with LQC observables 

To show it let us focus our attention to just one limit, 
say in the distant past. Due to the symmetry of the sys- 
tem the reasoning is immediately applicable also to the 
distant future one. Also, since the relation is of physical 
interest only when a± stay finite we restrict our studies 
to the states ^ which are localized in the weak sense, 
that is for which the expectation values and dispersions 
of the components of their WDW limit are finite (for 
finite </>). 

The forms (PTMl) . (jO^ of the LQC and WDW states 
and the relation between and its limit given by p.lOp 
and (|3.1ip imply immediately that 

11*^11 = 2||*||. (5.14) 

On the other hand, via the mapping ()2.15p these physical 
norms can be expressed as the appropriate kinematical 
norms on the surface ip = <j)o = const, which are given by 
(12. 5p and (|2.29l) respectively. This allows us to define the 
partial norms || • |1^ ^ ) as the restrictions of (|2.5p . (|2.29p 
to those points in the domains which satisfy ln|t;| > x 
for and In |u| < a; for '-' respectively. 

Consider now an arbitrary small e > and select the 
point Xo such that ||^!.^||(^ ^ j < Next define 

the function 

i{<j,)=Xo~{a/2){(j,-(l)o). (5.15) 

Using this function as a separator we can define the par- 
tial dispersions of the observables ln|{)|0 per analogy 
to the partial norms, that is restricting the domains of 
summation to the sets In \v\ > x{4)), In \v\ < x{(f>) respec- 
tively. In the similar way we define the dispersions 

of analogous observable acting on the state . They 
obviously sum up to the complete norms and dispersions 

ll*iP = (ll*IIJ)' + (11*11^)', (5.16a) 

= + (5-16b) 

m+lP = (ll*+IIJ)' + (ll*+ll^r, (5.16c) 

4 = (4)' + (^^)'' (5.16d) 

where (T0 is the dispersion of an observable (I2.26P and 
cr+ is defined via (I5.2bp . Using the known asymptotic 
behavior of the basis functions we can compare the limits 
of these dispersions as (j) —?' — oo. 

Let us start with cr^ . Applying the numerical estimate 
p.ip and taking into account, that the state ^' is localized 
in p^ we can write the asymptotic behavior of the wave 
function (for large negative (f>) 

*(t;,^) = Jf+(t;,^) + Oi(t;-5/2), (5.17) 

where the remnant Oi has the bound independent of (j). 
This, together with the fact that the value ^'^{v, (f>) de- 
pends only on In |w| — acf) and the convergence of the sum 



12 



over the relevant part of and integral over the domain 
ln|i;| > implies 

||*|i; = |||*+li; + 02(e-^(«), (5.18a) 
4=4 + O3(i(0)e-^('^)). (5.18b) 

On the other hand, the localization of in Injuj^ 
and the fact, that 

\x{(j))-x{(j))\ cx a;(<^), x{(j)) := \^^+f^ ' ^ ^ ^ 

implies that 

lim a^=0, lim ||*+||^i(0) = 0, (5.20) 

where the right-hand side equality is the consequence of 
the left-hand side one and ()5.19p . Furthermore, from that 
and (j5.16dp follows 

lim CT+ =(7+. (5.21) 

Combining together (l5TT4t . (I5.16cp . (I5.18bp and (j5.16al) . 
taking into account, that the part of "S/ contributing to 
the norms and dispersions is supported on ln|D| > 0, 
estimating the partial dispersion by the partial norm 

CT^< 11*11^ sup \\n\v\-xm (5.22) 

tie£°n[o,cxp(s(()i))] 

and applying (I5.19P we get 

lim ct7 = 0. (5.23) 

This resuh, together with (I5.16bp . (|5.18bp . ((5?T5|) and 
(|5.2ip finally implies 

lim CT0 — <^+- (5-24) 

Due to symmetry of the system, repeating the above rea- 
soning for the limit (/>—>■ cx) and the WDW component 
^ we obtain the convergence in the asymptotic future 

lim = cr_, (5.25) 

where (T_ is defined via (j5.2bp . Thus, provided the con- 
sidered LQC state ^ has finite dispersion inp^, the trian- 
gle inequalities (|5.10l) , (|5.12l) apply also to the asymptotic 
past and future limits of the dispersions of LQC states. 

VI. CONCLUSIONS AND OUTLOOK 

In this article we introduced the interpretation of an 
evolution of a universe described via Loop Quantum Cos- 
mology as the scattering process of the geometrodynam- 
ical (Wheeler DeWitt) one. Using this picture we ana- 
lyzed certain properties of the bounce in the model of 



a flat isotropic universe with the massless scalar field as 
well as in the isotropic sector of its generalization to the 
homogeneous but anisotropic spacetimes (Bianchi I). In 
these cases the LQC evolution is a process of transition 
of an ever-contracting (incoming) WDW universe into 
an ever-expanding (outgoing) one per analogy with the 
Klein-Gordon wave packet coming from the infinity, in- 
teracting with the nontrivial potential and being refiected 
in a scattering process back to infinity. The process is de- 
scribed by a scattering matrix, which explicit form was 
found and of which properties were investigated in detail. 
These properties were used to compare the dispersions 
of the observables In - scaled logarithmic volume at 
given "moment" of a scalar field (an internal time) of 
the incoming and outgoing states. It was proved that 
the dispersions satisfy certain triangle inequalities (j5.10[ 
15.121) involving also the spread of the logarithmic value 
of a scalar field momentum ln(p0/(?i\/G)). The derived 
inequalities are: (i) exact and, {ii) general, as they hold 
for any physical state admitted by the model. 

These results were immediately extended to the infinite 
past and future limits of the dispersions of true observ- 
ables acting on the genuine LQC states, as it was shown, 
that these dispersions converge in the asymptotic past 
and future to the appropriate dispersions of the, respec- 
tively, incoming and outgoing WDW states considered 
in the scattering picture. This convergence happens for 
every state on which the dispersion of is finite. 

The result reported above immediately implies that 
once the incoming state is semiclassical (that is it is 
sharply peaked in Injuj^ and \n{p^ / {h^/G))) so is the 
outgoing one and vice versa. In consequence any physi- 
cal universe of finite semiclassical in the infinite past 
is also always semiclassical in the infinite future. This 
conclusion is a precise confirmation as well a generaliza- 
tion of the results of [9] as it (i) holds for any physical 
state satisfying just a very reasonable condition of the 
finiteness of cTp^ , {ii) is valid also for the prescriptions dif- 
ferent that sLQC and which in particular are not exactly 
solvable, [Hi) extends immediately also to the isotropic 
sector of the Bianchi I model quantized in f20l] . 

In addition to confirming the robustness of the recall 
picture for the simplest models, the presented work con- 
stitutes the development of a new methodology, which 
application is not restricted to exactly solvable mod- 
els. Our method uses only the structure of the physical 
Hilbert space and asymptotic properties of its basis, none 
of which have to be controlled analytically. The only 
relevant requirement for application is sufficiently good 
control over states in the geometrodynamical analog of 
the LQC model under study. This allows in particular 
to easy extend our studies to more complicated isotropic 
models admitting well defined WDW limit, for example 
the ones with positive cosmological constant. For such 
models by construction our method seems to be better 
suitable to investigate asymptotic behavior (very large 
size) of the universe than the ones using the Hamburger 
decomposition [63| , as in this regime they reach the limit 
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of their applicability [63} . 

Presented methodology can be extended also to any 
system, in which one, while not having a good control 
over the large v limit of the LQC state, one still can ver- 
ify the semiclassicality preservation of its WDW analog. 
There, instead of comparing the incoming and outgoing 
components of the WDW limit at w — oo, one can 
compare them in the limit u — >■ 0, where they usually 
converge to the plane waves characteristic for the model 
we studied in this article. Such comparison alone is not 
sufficient to give any useful information about actual dis- 
persions in context of LQC, however it provides a con- 
trollable bridge between the WDW components. Then, 
once we are able to control (bound) the growth of dis- 
persions along the evolution of the WDW states, we can 
use this bridge to relate the dispersions at large v of 
and . Such method seems to be viable for example 
in the case of the universe with massive scalar field (i.e. 
the inflaton [13 )■ There, the WDW model is much easier 
to handle as, in opposition to the LQC one, it admits a 
good internal clock, evolution along which is generated 
by a self-adjoint operator. 

In certain, less precise sense an extension can be made 
also in the case of the recollapsing models. There how- 
ever, as in the physically interesting domain (universe 
sizes not bigger than the size at the recoUapse) the WDW 
limit of the state only approximately approaches the gen- 
uine LQC one and the real convergence happens already 
for the "tails" of the wave function, the parameters (like 
dispersions) of the WDW limits can resemble the analo- 
gous ones of LQC state only approximately, without ac- 
tual convergence. Therefore any found relation between 
the dispersions of the components of the limit can provide 
analogous relation for the LQC state only at the approx- 
imate level, up to some finite corrections. To get exact 
relations, like the triangle inequalities derived here, the 
precise estimates on those corrections need to be made. 
On the other hand, the presented scattering interpreta- 
tion is well fit to investigate in systems featuring quasi- 
periodic dynamics the phenomena like spontaneous co- 
herence in some epochs of the universe evolution. 

The methodology described in this article can be ap- 
plied not only to the isotropic models, but also to the less 
symmetric ones. Examples of such are Bianchi I models: 
vacuum ones or admitting the scalar field. In particular 
the analysis of the latter model, quantized via the meth- 
ods of [lal and [l3l will be presented in [63|. There, the 
semiclassicality preservation has a slightly weaker sense, 
as, due to nontrivial dispersion relation u){k) only the rel- 
ative logarithmic observables (A In jw^ |0)/(ln jwij^) have 
well defined finite limits as — ±{X). The results de- 
scribing the semiclassicality preservation in the vacuum 
case (quantized as in [31) are already presented in [To} . 
There, due to subtleties related with the choice of an 
emergent time the analogs of triangle inequalities con- 
tain additional terms, thus the relations are weaker than 
in the case with the scalar field. 
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Appendix A: Mathematical aspects of the analysis 

In this Appendix we discuss certain aspects of the anal- 
ysis, which, while are applied in our studies, are not di- 
rectly related to the topic of investigation. We also recall 
some general features of quantum mechanics which were 
used in the main body of the paper. 

1. The existence of self-adjoint extensions of 
Bianchi I evolution operator 

The evolution operator 0bi appearing in (I2.19p is a 
symmetric operator in the appropriate domain dense in 
Hkin- In order to define the physical Hilbert space it 
has however to admit at least one self-adjoint extension. 
Here, applying the analysis of the deficiency spaces [gl] 
we show, that this is indeed the case. 

The deficiency spaces are the spaces of kinemati- 
cally normalizable solutions -0* to the equation 

Qbi^^^ = ±i'ip^ ■ (Al) 

The existence of the self-adjoint extensions depends on 
the dimensionality of U^: if dim(?7"'") = dim(J7~) the 
operator admits the desired extensions. In particular if 
both the spaces are trivial, the extension is unique. 

At present the complicated structure of 0bi makes 
finding the solutions to (|Aip very difficult, however to 
achieve the task at hand one just needs to demonstrate 
the equality of dimiY^. To show that we construct a 
1 — 1 correspondence between the elements of these two 
spaces. 

Suppose, that ip^ is the element of that is it sat- 
isfies the equation (Ai, A2, ujOBi — jljV'''') = for every 
basis element (Ai, \2,v\. Expanding this set of equations 
in terms of ■i/;^(Ai, A2, v), acting on it with complex con- 
jugate and recomposing again one can see immediately 
(by inspection of the Eqs. (3.35)-(3.37) of 20]), that 7/^+ 
is the solution to (Ai, A2,w|0'^ -I- = 0, that is it 

belongs to . Since this reasoning can be also repeated 
in the opposite direction, the transformation 

R:U+^U-: [i?i/;+](Ai, A2, 1;) = 7^+(Ai, A2, 7;) (A2) 
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is a bijection. In consequence the dimensions of are 
indeed equal. 



the remnant parts will always sum up to a finite quantity 
we get 



2. Relation of norms of the LQC states and their 
WDW limits 

In most LQC quantization prescriptions considered in 
the literature the basis functions - normalized eigen- 
functions of G, can be evaluated only numerically. This 
is done by solving the difference equation p.22p for some 
chosen initial data ek{v = e). One does not know how- 
ever, for which value of efe(e) the solution is precisely 
normalized. Therefore in numerical studies one usually 
calculates the eigenfunctions which are not normalized, 
evaluates their norm and rescales them appropriately. 
However in many models, like the one considered here, 
the basis functions are normalized to Dirac S, thus the 
norm cannot be computed by purely numerical methods. 
Fortunately, the self-adjointness of 6 implies quite sim- 
ple relation between the norm of any eigenfunction of 
and the norm of its WDW limit. Since for given eigen- 
function this limit can be calculated numerically (see 
and Appendix |B]) , that relation provides sufficient data 
to normalize the LQC basis functions. Such method was 
implemented for example in '3, "7] and \TS\ . Although the 
discussed relation was applied already in those papers, 
due to lack of space its derivation was never presented. 
We show it here, since that relation is a key ingredient 
applied in the studies of Sec. IV Bl 

The derivation is essentially an estimate of the prod- 
uct (efc'|e/c) by the products (e_|_j,, |e_|_j.) via use of the 
asymptotic relation p.l[) . For simplicity we restrict the 
derivation just to the case e = 0, restricting the support 
of the eigenfunction to £q , although it immediately ex- 
tends to the remaining superselection sectors. The only 
difference in these sectors is the need to take into account 
both the limits v — >■ ±oo in some prescriptions. The in- 
termediate relations presented above are well defined in 
the distributional sense. 

Let us start with the orthonormality condition for the 
LQC basis functions 

X{k, k') := {eM = Y,'^^ek{v) = 6{k - fc'). (A3) 



She sum in the above equality can be spit as follows 

Xik,k') = Fi{k,k')+J2^^ekiv) (A4) 



where £i := >C(j^|t,>i and 



£+n[o,i] 



(A5) 



vanishes in the sector under consideration and is well de- 
fined function in the general situation. Applying the limit 
(j3.ip and taking into account, that the terms containing 



X{k,k') = F2{k,k')+r{k)r{k')x 

^^;^^^^(^) (A6) 

Ci 



E 



X > e" 

s,s'=±l 



The sum over £i can be now estimated via integral dv 
[66l |. The form (|2.31|) of e_|_j. implies that the correction 
due to this estimate is again well defined function of fc, fc', 
thus 



(A7) 



X{k, k') = Fsik, k') + -r(fc)r(fc')x 

/OO 

Knowing the form of e^^. and the relation [67| 



dxe{x)e^^-^\Uk)-^] (A8) 



we can evaluate the integrals in (IA7p . The result is 



X{k,k') = Fi{k,k') + -r{k)r{k')x 
8 

s,s' = ±l 



(A9) 



where F^ is again a well defined function. 

Taking into account, that k,k' e R+, thus the test 
functions integrated with the distribution have support 
only at positive k' we get the relation 



6{k-k') =X{k,k') ^Fiik.k') 



+ - cos(a(fc) - a{k'))r{k)r{k')5{k - k') , 

which can satisfied in the distributional sense only if 
r(fc) = 2. In consequence, given a WDW limit ^jj^ of the 
LQC basis function as defined in ()3.1|) . the following 
holds 



llV'fell =2V2||efc||. 



(All) 



3. Applicability of the stationary phase method 



In Sec. lV Al the stationary phase method was applied to 
approximate the integral p.20p determining the propor- 
tionality factor N . However, as the integrated function 
is singular in = and the integration is performed over 
the real line, the question whether the method can be 
applied there, is nontrivial. Here we show, that the con- 
tribution from the neighborhood of the singularity can be 
neglected and the method selected can be in fact applied. 
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Consider now the integral in (I3.20[) . It can be split 
onto two parts. 



L{k) := / dye^^'^'C"!^!-^) 
Jr 

= L+{k)+L-{k) + L{k), 
L^ik) / dyay)e''"^'''^«^-'"\ 
m / dy(l-^(2/))e-*'=('"l^l-^), 



(A12a) 

(A12b) 
(Al2c) 



where ^(y) is a smooth function equal to 

'O, \y-l\<P 

<( *, p<\y-l\<2p 

1, \y-M>'2p 



(A13) 



for some chosen small p. Stationary phase method can 
be safely applied to the term L{k). On the other hand 
one can show, that the remaining terms are bounded by 
the function 0{k^^). Indeed, integrating L^{k) by parts 
we get [eil 



k{l-e) 



lye 



£ = 

ik(ln(y)-y) 



Ml-.) 

1 d f y 
'"y-k^^y^d-yKT^y 



^-ik{\n{y)-y) 



where the integral in the 3rd line (denoted further as 
h^k)) can be rewritten as 



h{k) 



Hi-yy 



-ik{\n{y)-y) 



(A15) 



The form of both the integrals above and the function 
£,{y) allows us to estimate L^{k) as 



|i+(fc)|<| + | 



(A16) 



where are some constants common for all k. In conse- 
quence we obtain an estimate for the integral which is of 
the order 0(l/fc). The same technique can be applied to 
L~(k), giving the bound of the same order. 

The similar estimate can be derived for the integral 
(I3.23P determining N' by introducing the splitting anal- 
ogous to (jA12p . Then the term L'^{k) -an analog of 



L~^{k)- can be reexpressed as 

L'+{k)^ / dy^(2/)ln(y)e-''=(i-(^)-^) 



^ j^tf)£]^g-^fe(l„(e)-e) 
J ■ -1 

AyU'{y)'^e-^^^-^^)-y\ 

k l-y 

(A17) 

where the first righthand side integral (denoted as /i(fc)) 
equals 



i 1 

-e 



kl-y 

oo 



z 1 - y + ln(y) 
J? 

ik{ln{y)~y) (A18) 



^2/e(y)^(^e-'=('"('')-^). 

Applying again the integration by parts to the first inte- 
gral, we obtain the following estimate 



l^i(fc)l<| + |- 
As a result, we obtain an estimate on L'+(/c) 

l^'^WI<| + l + l> 



k 



(A19) 



(A20) 



where, as before, one can choose a common values of 
for all k. Analogously we arrive to the similar estimate 
on L'~{k) - an analog of L~{k). 



4. Relations between the dispersions and the 
correlations 

Here we briefly remind the well known relation between 
the correlation of two operators and their dispersions. 
That relation is general and holds for any quantum me- 
chanical system. While being quite basic, it seems to be 
often forgotten in the present literature in LQC, which 
sometimes may lead to an impression, that the correla- 
tions can grow uncontrollably. That relation is applied 
here in the derivation of the bound on the dispersion 
of the sum of the operators, which is in turn used in 
SeclVH 

Consider two operators A, B essentially self-adjoint in 
some dense domain in a Hilbert space T-L. Denote the ex- 
pectation values of these operators evaluated on the state 
^' S "H respectively as a and h. Denote also their disper- 
sions as a A and gb- The correlation E{A, B) between A 
and B is defined as the expectation value 



E[A, B) ((i - al)(B - 61)) + {{B 



bt){A-at)). 

(A21) 
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Applying the Schwartz inequahty to the states |$) := 
{A - al)|*) and \x) ■= {B - 61)1*) we immediately get 
the bound 



i\EiA,B)\^mx)\<\\m\x\\^<JA<JB. 



(A22) 



The use of Schwartz inequalities allows to derive analo- 
gous relations between at least some of the higher order 
components of the Hamburger decomposition (defined for 
example in [63|V Thus, it might provide a useful tool for 
the qualitative analysis of the physical evolution of the 
states described in terms of those momenta, in particular 
allowing for some control over their dispersion growth. 

Let us now consider the sum of operators A and B. Its 
dispersion on the state ^E* by definition equals 



^l+B = ^l + E{A,B)+al. 



(A23) 



Applying the bound (IA22[) and taking the square root of 
the above expression we arrive to the inequality 



(A24) 



which holds for any pair of self-adjoint operators and 
any physical state, independently of whether A and B 
commute. 



Appendix B: Numerical aspects of the analysis 

In this appendix we describe the actual numerical tech- 
niques used to investigate the properties of the phase 
shifts a{k) defined in p.l|) . which were presented in 
Sect. Hire] and applied in Sect. FOl 

The eigenfunctions ipk — f{k)ek can be calculated di- 
rectly from (|2.22l) via iterative methods described in Q- 
Since (|2.22p is of the 2nd order, to determine the eigen- 
function uniquely one in principle needs to provide its 
initial values at two consecutive points of Cg. However 
for the prescriptions and the superselection sectors con- 
sidered in this paper the value — 8) is determined 
uniquely by il)k{v = 4) (see 0,11, [HI for the details), thus 
the latter one is the single initial datum. A main disad- 
vantage of this method is, that at the level of providing 
the initial data it is not possible to normalize the func- 
tion. It can be done only via application of the relation 
(IA11[) . once V'fc is calculated on some large domain of v 
and its WDW limit is found. Fortunately the phase shifts 
a defined in (jSTj are not sensitive to the normalization, 
thus in our studies '0fe(4) (or equivalently f{k)) can be 
set arbitrarily. We fix it choosing 



MV^4:) = 1. 



(Bl) 



Once f/jfe is evaluated in some domain jCq fl [0,Um] 
(where vm ^ |^|), its WDW limit can be found for ex- 
ample via an extrapolation to t; — > (X) of the coefficients 
defined in p.Sp . In actual simulations to find a a more 
stable and faster converging method was implemented. 



Namely the vectors Xk {v) were decomposed analogously 
to (prO)) . that is 



Xk{v) = r{k,v) ( g_i„(fe,„) 



(B2) 



and the components Q;(fc, v) were calculated only at 
the points u„ = 4n, where the signs sgn(7/)fe(D„)) ^ 
sgn('0fc (■(;„ — 4)) with use of the observation, that the 
function 



(B3) 



il^iv) := r{k,Vn) 

X [e'^^^e^'^-^e^w) -He"*'="('=^^")e_fc(?;)] 

vanishes only if 

fcln(i;) + a(A;,u) = 7r/2 + TO7r, m e Z. (B4) 

The position of the roots (denoted further as Vo^n — 
exp(a;o,n)) was identified via the linear approximation in 
X = ln(u) 



ln(wo,n) = ln('y„ - 4) 



-In 



(B5) 



ilkiVn) - 1pk{Vn - 4) ' 



Once the set of Q!(fc,Wo „) was determined, the limit 
a{k) was found via a polynomial extrapolation (Neville 
method) with respect to the variable y = exp(— Xq) at 
2/-0. 

Such method has been applied for example in 7] to 
evaluate a as a basis for the analytic approximation fur- 
ther used for construction of the initial data correspond- 
ing to the symmetric states. However, while it allows to 
calculate a quite precisely, it cannot be directly applied 
to find its derivatives, as it can be differentiated only 
numerically, whence 

(i) the numerical errors accumulating over the steps in 
the computation of a would not allow to derive suf- 
ficiently precise differential limit for small displace- 
ment of k, and 

(ii) for larger displacements of k one could not exclude 
the existence of small "large frequency" terms which 
in principle, while giving negligible contribution to 
a itself, could give a considerable one to its deriva- 
tives. 

Here, to avoid the above problems, we use a slight mod- 
ification of the method under discussion, namely we find 
a' by differentiating (IB4[) over k 

(B6) 



[dka\{k,v) ^ -{xo + kx'o), Xo^\n{vo), 



where the values x'^ „ of the derivative x'^ are in turn 
determined via differentiating (|B5P 



I'n -4 
^fc(^n)V'fc('i^' 



(B7) 



(V'fc(Wn) - ^k{Vn - 4))^ 
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The derivatives 4''i.{vo,n) appearing in the above expres- 
sion can be in turn evaluated via an iterative equation 
resulting from a differentiation of ()2.22|) . that is 

^'kiv + 4) = -^[(w' + fo{v))^'k{v) + ^gg^ 
-UvW,{v-A)]. 

Up to the inhomogeneous term proportional to ipkiv) 
that equation is a complete analog of (j2.22p . thus the 
initial value problem is equivalent as well, once ipk itself 
is known. Therefore the initial condition (IBip specifies 
uniquely not only 7/;^ but also V-'j.- 

The equations jHH), dUl, dlZl and dBSj allow one to 
determine [9fca](fc, Wo,n) corresponding to finite u„. Then 
the limit [dkO:]{k) is found via the polynomial extrapola- 
tion in the exactly same way as a{k). 

To find [d1a]{k) one proceeds in the same way, first de- 
termining [dla]{k,Vo,n) analogously to [dka]{k,Vo,n) via 
the system of equation build by differentiation of (IBip . 
(|B6|) . (|B7)) and (|B8| and computing the final limit by the 
polynomial extrapolation. 

In the actual simulations performed the sequences 
{xo,n}i of the roots were searched for in such a way, that 



an zth root in the sequence is the maximal root within 
the range [0,2~^vm]- The bound vm ranged from 10^ to 
2 • 10^. Depending on the number of the roots actually 
found the order of the polynomial used to find the final 
limits varied from to 2. The range of k investigated 
numerically was [10~^, 10®]. 

The results are presented in Figs. [T] to [3] As we can 
see in Fig. [1] the chosen methods allow to calculate dkOi 
quite precisely within the entire investigated domain. Its 
precise form varies for small k depending on the par- 
ticular prescription, whereas for large k the derivative 
quickly approaches — ln(/c). The rate of approach p.26p 
is confirmed for k < 10^ (see Fig. [5^), however due to 
the numerical errors the results are no longer reliable for 
larger k. Nonetheless the results converge also there as 
Vm increases, which is shown in Fig. [^b- 

The 2nd order derivative could be evaluated reliably 
only for k < 10, as it is more sensitive to the numerical 
errors. The result is presented in Fig. [5^ showing the 
rescaled derivative kd^a. From there we observe that it 
quickly approaches its limiting value kd^a — >■ — 1 deduced 
from Fig. [T]i, never exceeding it [6^. As in the case of 
the remnant of dkCt the numerical results converge as vm 
increases, which is presented in Fig. [SJj. 
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FIG. 1: The behavior of dkce for different quantization prescriptions compared against its large k hmit a'{k) = — In(fc). Figure 
(a) shows the shape of the derivative in the entire domain investigated numerically, whereas (6) presents the zoom for small 
values of k. 
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FIG. 2: The values of the rescaled remnant function /?(fc) :— (a'(fc) + ln(fc))fc^ calculated numerically within the domain of 
reliability of the applied method. Figure (a) shows the result for different prescriptions whereas (6) presents the dependence of 
calculated /3 on the choice of vm on the example of the APS prescription. One can see that as vm increases /3{k) converges. 




FIG. 3: (a) The rescaled derivative kd^a for different quantization prescriptions shown in the domain, where it could be 
determined numerically in the reliable manner, (fo) The values of kd^a for APS prescription evaluated on different sets of 
probing points with the upper bound ranging from 10^ to 2 • 10*. The numerical solutions converge as the bound increases. 
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